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The study of the broadband bosonic communication 
channel has been one of the first applications of quan- 
tum communication theory |lj, |i2,] . The basic result of this 
effort has been the determination of the ultimate limits 
posed by quantum mechanics to the rate at which classi- 
cal information can be reliably transmitted through the 
channel in the noiseless case. In this context, the classical 
capacity C was shown to be proportional to the square 
root of the input power. Here we generalize these results 
by extending the analysis to noisy configurations and to 
other channel capacities such as the quantum capacity Q 
(the amount of quantum information that can be reliably 
sent through the channel) or the entanglement assisted 
capacity Ce (the amount of classical information that can 
be reliably sent through the channel in the presence of 
an infinite amount of prior entanglement between sender 
and receiver). In particular, we study various types of 
noise: loss (where there is a probability 1 — that a pho- 
ton is lost in the transmission line), loss with white noise 
or thermal reservoir coupling, and a "dephasing" channel 
(in which the average number of photons is preserved, but 
some phase correlations are lost in the transmission). In 
this context we determine the value of Ce as a function 
of the input power and show that the square root de- 
pendence applies also to most of these channels. For the 
other capacities we provide some bounds that establish 
the same dependence. A sketch of the results obtained 
is summarized in table even though implicit equations 
for all the capacities (or their lower bounds) have been 
obtained, in most of the cases numerical methods have 
been employed to derive their values as a function of the 
channel parameters. 



We start by introducing the model of the channel and 
of its noise sources in Sec. 13 We introduce capacities 
Ce, C and Q and the Lagrange procedure that is used to 
evaluate them for the broadband channel in Sec. |nj The 
remaining sections are devoted to the analysis of the lossy 
channel (Sec. IIII|I . the white noise channel (Sec. lIV|l . the 
thermal noise channel (Sec.EJ and the dephasing channel 
(Sec. ED. 
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TABLE I: Capacities calculated in the paper for different 
noise models. 

I. GAUSSIAN BOSONIC CHANNEL 

The prototype of a high capacity communication chan- 
nel is an optical fiber, where time or frequency multiplex- 
ing (or hybrid strategies) are used to send information. 
From a fundamental point of view, such a communica- 
tion line is described as a broadband bosonic channel 0] ■ 
In the present paper we analyze the performance of this 
channel in the realistic scenario of non-perfect transmis- 
sivity, i.e. the possibility that photons can be lost during 
the communication or that they can be replaced by pho- 
tons coming from external noise sources. The analysis 
is complicated by the fact that, for some capacities, it 
is not known whether the additivity property holds, i.e. 
whether entangling successive uses of a noisy channel may 
increase its transmission rate 3j. 

Without loss of generality, we will assume that for 
each frequency in the channel only one polarization is 
used to transmit information, i.e. no frequency degener- 
acy is present. The quantum description of this channel 
is obtained by coupling each mode to a noise reservoir 
with beam splitters that have transmissivity equal to the 
quantum efficiency rjj of the jth mode, i.e. 

= v% + V^~-Vj bj , (1) 

where aj, a'j and bj are the annihilation operators of the 
input, output and noise modes respectively. The loss 
map J\fj for the jth mode arises by tracing away the noise 
mode bj and the global loss map A/" is the tensor product 
{^jA/i. Notice that for rjj = 1 the noise reservoir is 
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decoupled from the transmission line: this describes a 
noiseless channel where N — Different types of noise 
can be described depending on the initial state of the 
reservoir modes hj. We will analyze the case in which 
the reservoir is in a separable Gaussian state of the form 

nib) _ 
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where z is a real bidimensional line vector, Aqj = qj — 
(qj), Apj = Pj — {pj), with qj and pj the quadratures 
q^ = ,JKj2 [b] + bj) and Pj = i^fhp. (&] -bj). For the 
situations in which we are interested, {qj) ~ {pj) = and 
the correlation matrix Bj in Eq. has the form 



B,=h 



Nj + 1/2 






Nj + 1/2 



(3) 



where Nj is the average number of photons in the noise 
mode bj. With this choice of pib\ the CP-map Af de- 
scribes a Gaussian channel, i.e. it transforms Gaus- 
sian input states (in the modes Uj ) into Gaussian output 
states (in a'j). 

Four different noise models will be analyzed in detail 
in this paper. The simplest one is a purely lossy channel 
in which the photons in the jth mode have a probability 
1 — rjj to be lost during transmission. It is described 
by taking Nj — for all j, i.e. by taking pib) in the 
vacuum. For optical communications this is the most 
interesting situation, since thermal photons are negligible 
at room temperatures. Another interesting case is given 
by choosing Nj = N for all j, which describes an added 
white noise to the transmission. On the other hand, by 
choosing 



1 



ghLOj/{KT) 



(4) 



(with LOj the jth mode frequency) we can describe the 
effect of coupling the communication line to a thermal 
reservoir at temperature T. Of course, in the limits N — )■ 
or T — > both the white noise and the thermal channel 
reproduce the lossy channel. The common trait among 
these three noise models is the fact that they can be 
parametrized as 



Nj=Nv{cOj/u,) 



(5) 



where N describes a characteristic number of photons 
in the transmission and uj describes (through an appro- 
priate function v) an eventual characteristic frequency 
of the channel. The parametrization ^ will be useful in 
deriving some scaling properties that simplify the deriva- 
tion. A final noise model we will analyze is a non-linear 
noise mechanism where the average photon number of the 
reservoir Nj is a function of the average photon number 
in the message. This type of model is well suited to de- 
scribe situations in which the noise is due to the action 



of some active third party (e.g. an eavesdropper) who is 
tampering with the transmission. In particular we will 
analyze a sort of dephasing channel where the average 
photon number in each transmission mode is preserved, 
even though some phase correlation is lost. Because of 
the non-linearity of this noise, the parametrization Q 
does not apply to the dephasing channel, but most of the 
general formalism developed for the other models can still 
be used in this case. 



II. CAPACITIES 

In this section we introduce the three different channel 
capacities that will be analyzed in the paper. 

Entanglement assisted capacity. — The one channel 
capacity that is known to be additive 0, 0] even in the 
presence of noise is the entanglement assisted capacity 
Ce- It is defined as the number of bits that can be re- 
liably transmitted per channel use in the presence of an 
unlimited quantity of prior entanglement shared among 
the sender and the receiver. This quantity gives a sim- 
ple upper bound to all the other channel capacities and 
is conjectured to provide an equivalence class for chan- 
nels Analogously, one can define the entanglement 
assisted quantum capacity Qe that measures the num- 
ber of qubits that can be reliably transmitted per channel 
use in the presence of an unlimited quantity of prior en- 
tanglement. Using teleportation and superdense coding, 
it is easy to show that Qe — Ce/2, so that only one of 
these two quantities needs to be determined 0]. 

Taking advantage of its additivity property, the entan- 
glement assisted capacity of a multimode channel can be 
calculated as 



ued capa 
Ce = max I I{J\ft, ft) \ , 



(6) 



where "Hj is the Hilbert space of the jth mode in the 
channel and I{Mj , Qj) is the quantum mutual information 
defined as 8] 

/(AG, ft) ^ ^(ft) + ^(AG[ft]) - S{{N, ® l)[<i>,J) , (7) 

(with S{gj) — — Tr[ft logj ft] the Von Neumann entropy 
and a purification of the mode input density matrix 
ft). The maximization © will be performed only using 
the states ft that satisfy the average energy constraint 



(8) 



where ujj is the frequency of the jth mode and A^ = 
Tr[a|ajft] is its average number of photons. This 
constraint is fundamental: without any restriction the 
bosonic channel would have infinite capacity since the 
Hilbert space that the sender could use for encoding 
would be infinite-dimensional. The energy constraint in- 
troduces an effective cut-off in the dimension of the cod- 
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ing space j^. This, of course, mirrors any realistic situa- 
tion in which the energy available for the transmission is 
always finite. 

Since we are dealing with a Gaussian channel M we 
can apply the Holevo- Werner theorem which states 
that I(Afj,gj) reaches its maximum over Gaussian in- 
puts. Moreover, for the noise models we analyze, squeez- 
ing the input to the jth mode does not increase its quan- 
tum mutual information if the energy of the mode is fixed 
(see App. 0for details). Hence, the maximum value of 
I(Nj,Qj) is given by an expression csiNj, Nj,rij) that 
depends only on the number of photons Nj , on the noise 
parameter Nj, and on the quantum efficiency rjj. The 
explicit form of cb is evaluated in Add. 1X1 and is 



(9) 



D,+N,-m-l\ rDj-Nj+N'-l 



where 



(10) 



is the average photon munbcr in the jth mode at the 
channel output and 
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1) 



(11) 



{x + 1) log2(a; + 1) — a;log2(a;) for a; 7^ 

for a; = 







In terms of c^, Eq. lO becomes 

Ce = niax I ce{Ni,Ni, rg)^ , 



(12) 



(13) 



where the maximization must be performed on the iVj-'s 
that satisfy ©. In Sec. Ill Al wc will calculate explicitly 
the right hand side of H13I) . First, however, it is conve- 
nient to introduce the other channel capacities in order 
to underline some common features. 

Classical capacity. — The classical capacity C mea- 
sures the quantity of bits that can be sent reliably 
through the channel per channel use (without assistance 
of prior entanglement). For the noiseless broadband 
bosonic channel {rjj = 1) it has been shown P, |^ that, 
under the energy constraint ©, 



C 



T Rc 



(14) 



where T is the transmission time, H. = ®iHi is the 
Hilbert space of the multimode channel {Hj being the 
space of the jth mode), and 




(15) 



is the classical communication rate in terms of the input 
power V = £/T. 



In the presence of noise, however, a recipe to calcu- 
late C involving only single uses of the channel, as in the 
case of Eqs. and H14|l. is not known. It could be that 
entangling successive uses of the channel the amount of 
information transmitted is increased 10, llj. This would 
require to consider input states in the Hilbert space 7i®" 
pertaining to n successive uses of the channel. The es- 
timation of C is, hence, a daunting task. Nevertheless, 
a simple lower bound for it can be obtained consider- 
ing unentangled coding procedures, where the sender is 
not allowed to employ codewords which entangle differ- 
ent channel uses. In the multimode channel, a further 
simplification consists in considering coding procedures 
where entanglement among the different signal modes aj 
is forbidden. This yields the inequality 



C ^ max 
Pj(Ai),o-j(/i; 



{^AV. (p^(A*),'T,(/i))} , (16) 



where Qj = J d/i pj(fi)aj{ii) describes a message in which 
the /ith "letter" encoded in the density operator aj (/i) G 
Hj has probability density Pj(^) to be sent in the jth 
mode and where 

AV, ^ S{M,[g,]) -Jd^i p,{^,)SiM,[a,i^,)]) (17) 

is the Holevo information. Unlike the case of Ce dis- 
cussed in the previous section, it is not known whether 
the maximum of Eq. (|16|l can be evaluated working only 
with Gaussian states. However, adopting this strategy 
one still obtains a tight lower bound for C 0. Thus, 
we evaluate Xj\f- {pj (/i), Uj (/i)) for the jth mode using co- 
herent states crj{fi) — weighted with Gaussian 
probability distribution Pj{fi) — exp[~\fi\'^ /Nj] / {nNj) . 
Selecting this encoding we are assuming that, as in the 
case of Ce, squeezing does not increase the unassisted 
capacity if there is an average energy constraint on the 
input state 

nil El El- With this choice, Eq. ^ can 
be written in a form analogous to H13|) . i.e. 



C ^ max|^fc(Ar, 



N. 



(18) 



where the maximum must again be taken under the av- 
erage energy constraint (jHJ . The function k is calculated 
in App. IXIand is given by 



kiN,,N„7j,) = giN'^) - g{{l - f^,)N,) 



(19) 



with iVj defined in Eq. (P|l . 

Equation H18|l establishes a lower bound for the classi- 
cal capacity C. Simple upper bounds for C are given by 
the entanglement assisted capacity Ce of Eq. (|33|l and 
by the noiseless classical capacity T Rq of Eq. ((TUl . 

Quantum capacity. — The quantum capacity Q of a 
channel is the number of qubits that can be sent reliably 
through the channel per channel use. For the noiseless 
case one can show that Q = C, i.e. for bosonic channel 
with j]j = 1 one can show that Q = T Rq. As for the 



4 



classical capacity, an expression involving only single uses 
of the channel is not known for noisy channels: again it 
could be that the entanglement of successive channel uses 
might increase Q [3 ■ Also here we will consider the 
lower bound obtained by excluding all the coding proce- 
dures that make use of entanglement among successive 
uses of the channel or among different modes. This pro- 
vides the inequality 



where 



Q ^ niax I ^ J(M, , 



J(AA, q) = /(AA, q) - S{q) 



(20) 



(21) 

is the coherent information [11 m. Equation ^ IS a 
consequence of the fact that random quantum codes can 
convey a number of qubits equal to the coherent informa- 
tion of the channel (if it is greater than zero) — for 
the particular case of the Gaussian channels this same re- 
sult was proved also in In evaluating the right side 
of Eq. (|20|l . we will employ Gaussian states: here the 
Holevo- Werner theorem does not apply and this choice 
will further lower the bound H2U|) on Q. Moreover, we 
will restrict the analysis to non-squeezed inputs. These 
considerations allow us to write Eq. (|20|l as (see App. 



Q ^ niax^(7(iVi,iVi,77i) 

* 7 



(22) 



where the maximization must be performed under the 
energy constraint © and (see App. ^ 



D, + N, -N'-l 



(23) 



with N'j and Dj defined in (|10|l and respectively. 

An alternative lower bound for the quantum capacity 
Q can be obtained by observing that the definitions of 
Ce in © and of / in (Tjl imply 



Ce — max 

gen 



[j{M,q) + S{q)^ s$ mjix{j(AA,e)} 



+ max {5(e)} ^ Q + max{5(e)} , 



(24) 



T Rc 



which for the broadband channel gives Q ^ Ce 
by employing Eq. I|14() . 

Equations H22|l and (|24|l give two lower bounds to Q. A 
simple upper bound is given by the entanglement assisted 
quantum capacity Qe = Ce/'^- 



A. Lagrange multiplier procedure 

In order to determine the values of Ce and the lower 
bounds for C and Q given by Eqs. and (|^ . 



one needs to perform maximizations of the form 



W — max / 



{j2w{N,,N,,Tj,)} , (25) 



under the constraint given by Eq. ((H)). In Eq. if^ . the 

quantity W represents Ce or the lower bounds for C 
or Q depending on whether w is equal to ce, k or q 
respectively. The Lagrange multiplier procedure is well 
suited to perform these constrained maximizations. It 
amounts to finding the values of {Nj} which solve the 
equations 



d 



1 



E 



LO,N,, 



0, (26) 



where l/(illn2) is the Lagrange multiplier that must be 
chosen to satisfy the constraint JSJl after having solved 
Eq. (EH)). The explicit expressions of Eq. (|26|l for the 
three capacities are reported in table |nj These equations 
are in general difficult to solve. A first useful simplifi- 
cation is to assume that all the modes have the same 
quantum efficiency, i.e. rjj — rj for all j. Even though 
this is a strong assumption, it is still a good description 
for broadband channels that have a wide spectral trans- 
mission window. Under this approximation, it is easy 
to verify that Eq. H2t)|) has solution that depends on the 
mode frequency ujj and on the noise parameters cu and 
N of Eq. © as 



n,=:f 



(27) 



To calculate il, the energy constraint ijHl can be written 
as 



duj ^ /uj to - 



2 /-oc- 



SU! Jo 



dx X J- 



{^^^'^'V) , (28) 



where the sum over the mode index i is approximated 
with an integral over the mode frequencies, under the 
assumption of small minimum frequency interval Sio of 
the channel. This last quantity determines the minimum 
time T = 2Tr /Suj needed to transmit a signal in the chan- 
nel. In order to solve H28|l in terms of il, it is useful to 
introduce the adimensional parameter i/q ~ uj/^l. Thus 
we find 



27rP 



1/2 



where V = £ jT is the average input power. 



/(2/o,iV,7;) 



dx X T{x,yo,N,r]) , 



(29) 



(30) 
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lower bound 


V 77,iV, + (l-77,)ivJ 




lower bound 





TABLE II: Lagrange equations deriving from 1)26^ for the different capacities. The functions A'^' and Dj are defined in pOjl and 
IIH respectively and Aj = [(1 — 3r]j)Nj + (1 — -qj) + (1 + rjj)N'j\/ Dj. Notice that the equation pertaining to C can always be 
solved analytically. 



and j/o is determined by solving (with respect to y) the 
equation 



(31) 



If the noise reservoir does not have a characteristic fre- 
quency u) (as in the case of the loss, white noise and 
dephasing), the derivation simplifies since neither T nor 
/ depend on the parameter y^: Eq. does not apply 
and is already determined by Eq. (|29|l . However, for 
the sake of generality, we can include also these last cases 
in the above formalism by assigning to them (D = and 
2/0 = 0. 

The value of W is finally obtained using Eqs. (|27|l and 
(I29|l to evaluate the sum 12511 . i.e. 



(32) 







duJ ( ^ ( ^ ^ ir^ \ I '-^\ 



where the parametrization Q was employed. Apart from 
corrections of order 1 /T (coming from the approximation 
of the mode sum with the frequency integral), Eqs. 
and H32|l imply that 

M^ = ri?c W(yo,iV,77) , (33) 

where Rc is the noiseless classical rate of Eq. (I15f) and 

In 2 \ 3 



W(2/o,^,r/) 



TT V 2/(yo,iV,77) 



(34) 



The quantity W is a proportionality factor that charac- 
terizes the dependence of W on the noise parameters ?7, 



IS! and cD. Even though it is in general difficult, if not im- 
possible, to analytically evaluate the expressions T and 
W, one can still provide numerical solutions for these two 
quantities as will be shown in the following subsections. 

When there is no characteristic frequency uj (as in the 
cases of loss, white noise and dephasing), then ?/o = and 
Eq. tells us that W depends on the input power V 
only through the classical communication rate Rc- This 
means that, in these cases, the capacities of the channel 
(or at least their bounds) are proportional to the square 
root of V just as the noiseless channel classical capacity 
of Eq. (|14|l . On the other hand, when a characteristic 
frequency w does exist (as in the case of the thermal 
noise) , then W depends on V also through the parameter 
j/o, which, according to Eq. (|31|l is a non-trivial function 
of up' jV ■ However, for fixed value of this ratio, the \fV 
proportionality still applies. 



III. LOSSY CHANNEL 

The simplest channel is the lossy channel where the 
reservoir is in the vacuum and iV = In this case, the 
Lagrange equation solutions T of (|27|l for the three ca- 
pacities are functions only of Wj/fi and r\. Thus, Eq. H33fl 
gives 



Ce = r Rc Civ) , 
C ^ T i?c ICiv) , 
Q ^ TRc Qiv) , 



(35) 
(36) 
(37) 



where the functions C, JC and Q take the place of W in 
Eq. (|34|l by replacing w with ce, k and q. These func- 
tions are plotted in Fig. ^ In all these three cases, the 
dependence on the input power is given by the VV term 
in Rc- The alternative lower bound for Q of Eq. H24|) 
becomes 



^TRc [C{r,) - 1] . 



(38) 
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Both this bound and the function Q{r]) of Eq. H37() are 
positive only for 77 > 1/2. This reflects the fact that for 
?7 ^ 1/2 the quantum capacity Q is nuU. A simple argu- 
ment based on the no-cloning theorem 18] is sufficient to 
prove this, as in the case of the erasure channel 19] . In 
fact, assume that Q is positive for 77 < 1/2, and suppose 
that a third party collects all the photons lost during 
the transmission: to him the channel would appear to 
have a quantum efficiency (1 — 77 > 1/2) greater than the 
one of the receiver. If Q > 0, both he and the receiver 
would be able to reconstruct the quantum information 
sent through the channel reliably, thus violating the no- 
cloning theorem. Interestingly this rj = 1/2 bound for 
the lossy channel has been observed also in tomographic 
quantum state reconstruction, where the effect of the loss 
can be deconvolved from the reconstruction only when 
the detection efhciency is bigger than 1/2 poj . 

Notice that, while C and Q must be computed through 
numerical methods (a partial analytic characterization of 
C{ri) is provided in Q, where, in particular, it is shown 
that C(l) = 2C(l/2) = 2), it is possible to determine an- 
alytically the value of /C. In fact, the Lagrange equation 
for C (see table lll|) has solution 



oWj7(o»)j) _ I 



(39) 



When all the rjj are equal, the value of O can be calcu- 
lated directly through the energy constraint {yo = 
since there is no characteristic frequency lu). In particu- 
lar, since the function / of Eq. (|30() is 



(40) 



we have Q. = [12V / [■nhr])]'^/'^ . To evaluate JC through 
Eq. H34|) . we need also the term 



f 

Jo 



ev -1 



7r^77 
3 In 2 



(41) 



Replacing gUJl and (glj in Eq. (PH), we finally find 

/C(77) = ^ . (42) 

Notice that, from Eq. lfn|) it follows that in Eq. 
equality must hold for 77 = 1 and the results of 
reobtained. 



IV. WHITE NOISE CHANNEL 

As in the case of the lossy channel, also the white noise 
channel has no characteristic frequency (D, in fact all noise 
modes contain the same average number of photons Nj = 
N. This means that the solutions T of the Lagrange 
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FIG. 1: Plot of the functions C(j]), ¥i(j]) and Q(j]) of 
Eqs. 1351 - 137II for the lossy channel. The function C(r]) char- 
acterizes the entanglement assisted capacity Ce (upper con- 
tinuous line). C(77), and the value 1 (that corresponds 
to the noiseless classical capacity 1141 ') bound C which is re- 
stricted to the hatched region. C{ri)/2 and 2(7?) bound the 
quantum capacity Q which is restricted to the cross-hatched 
region. The alternative lower bound for Q of Eq. (1381 is given 
by the dotted line, which is almost indistinguishable from 
2(77). The quantum capacity Q is null for 77 < 1/2 (dashed 
line). The fact that, for the lossy channel C(l) = 2 is a sig- 
nature of the superdense coding effect: prior entanglement 
allows to double the channel capacity [2H . 



equation are only functions of tOj/^l, N and 77. Thus 
Eq. |j3Sl) gives 

Ce = TRcC{N,7j) , (43) 

C ^ T RcJCiN,v) , (44) 

Q > TRc Q{N,ri) , (45) 

while Eq. ^ becomes Q ^ T Rc [C{N, 77) -1]. As in the 
previous case, Eqs. H43|I " H45() display the square root de- 
pendence of the capacities on the input power V through 
Re- The functions C, /C and Q for the white noise chan- 
nel are plotted in Fig. [3 for different values of N. Some 
examples of the numerical solutions !F{ijJi/VL, N, 77) of the 
Lagrange equations are plotted in Fig. 13 

As for the lossy channel, also here an analytical ex- 
pression for /C exists. In fact, the Lagrange equation for 
C (see table has solution 



N . (46) 



Since Nj represents the average photon number in the 
jth mode, the solution H46() can be used only when Nj ^ 
0. This condition is satisfied only if the frequency of 
the mode Uj is lower than the cut-off frequency LUmax — 
T] Q s, where 



In 



1 



(1 - Tj)N 



(47) 
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FIG. 2: a), b), c) Plot of the functions C{N,ri), K.[N ,7^) and 
Q{N, if) of Eqs. lH^- lUSt for the white noise channel with 
increasing along the direction of the arrows. The dotted lines 
represents the case iV = from Fig. d) Plot of the bounds 
of the classical (hatched region) and quantum (cross-hatched 
region) capacities for TV = 1. From top to bottom, the curves 
are C, K,, C/2, Q and the alternative lower bound C — 1, which 
in this case is practically coincident with Q. 




FIG. 3: Examples of the solution iV, = T{uj/n,N,ri) of the 
Lagrange equations for C'e (continuous line), for C (dotted 
line from Eq. 1461 '). and for Q (dashed line). 



(again we have assumed rjj = rj for all j). For ujj ^ 
i^max, Eq. cannot be used (it gives a negative Nj) 
so that we assume Nj = 0. This physically corresponds 
to not sending any photons in the high frequency modes: 
it would be too expensive in energetic terms to contrast 
the noise in these modes. With this choice, from Eq. H3U|) 
we obtain 



dx X 



1 



{l-r,)s^N/2\ , 



'-N 



(48) 



where 

nx)^j\yj-^. (49) 

To calculate K. we also need the integral 

dx k {T{x,N,r]),N,r)) 

dx [g {f^T + (1 - 7j)N) - g{{l ~ 7j)N)] 
1 



rjs 



V 



dx g 



V - 1 



sg{{l-7j)N) 



(50) 



Repjacing (gHl and (ISOJ into Eq. (PU the function 
JC{N , rf) is determined. Notice that in the limit iV — > 0, 
it is possible to show that the function K, converges to 
y/yy so that one reobtains the results of the lossy channel. 



V. THERMAL NOISE 

Let us now analyze the case of thermal noise. This 
noise model does have a characteristic frequency Co = 
KT/h which depends on the bath temperature T. Since 
(I) 7^ 0, we need to solve ?/o from Eq. H31(l which clearly 
implies that yo is a function of the ratio between the 
square of the temperature T and the power V. In this 
case, N = 1 and the expressions for the capacities are 



Ce = r Rc C{y'^^ ,r^) , 
C > TRcICiy^,v) , 
Q ^ TRc Qiy^,v) , 



(51) 
(52) 
(53) 



Ce „,c g^j-^^ yQ solutlous of Eq. (|31() for 



where yQ-^^yS 

the respective capacities. Moreover, the alternative lower 
bound of Eq. ^ gives Q ^ T Rc [C{y^'',r]) - 1]. The 
presence of the terms yo complicates the dependence of 
the capacities on the input power V. However, once the 
ratio H{KT)'^/V has been fixed, the usual dependence 
on the square root of the input power applies. Some 
numerical plots of C, /C and Q are given in Fig. 0]as a 
function of r/ and of the temperature T. Some examples 
of the corresponding Lagrange equation solutions J-' are 
given in Fig. [S] 

Again the function K. for the lower bound of the classi- 
cal capacity can be solved analytically. We will find that 
below a critical temperature the solutions Nj of the 
Lagrange equation, i.e. 



(1 



Vj)/Vj 



1 



(54) 



are valid for all frequencies LUj. On the contrary, when 
T > Tc a, cut-off frequency arises above which (as in the 
case of the white noise channel) it is convenient not to 
send photons. 
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FIG. 4: a), b), c) Plot of the functions C(yf^'^ ,r)), /C(yf ,7?) 
and Q{y§ ,1]) of Eqs. I|51^ - I|53|l for the thermal channel with 
temperature T increasing along the direction of the arrows. 
The dotted Unes represents the case T = from Fig. Q In 
the plot c), the low temperature regime T < Tc (continuous 
lines) is obtained from Eq. (1571 . while the high temperature 
regime (dashed lines) is obtained from Eq. 16UII . d) Plot of the 
bounds of the classical (hatched region) and quantum (cross- 
hatched region) capacities for Rt/Rc ~ -41. From top to 
bottom, the curves are C, K., C/2, Q and the alternative lower 
bound C — 1. 




FIG. 5: Examples of the solution Nj — T{lu/Q, N ,1]) of the 
Lagrange equations for Ce (continuous line), for C (dotted 
line from Eq. EHDi sjad for Q (dashed line). It is possible to 
show that the solutions of the Lagrange equation for Q have 
two cut-off frequencies for low and high lj. In this graph only 
the first one is evident. 



obtained from the energy constraint (|29(l is 
^ 6 In 2 / ' " ; 7 

Replacing Eqs. and (EHl) in Eq. (I23J), we obtain 



(56) 



where the dependence on y§ derives from (|55() and where 



A{x) = In 2 / dy g 



e 'y - 1 



(58) 



Equations H54|l and (|5(j|) are consistent (i.e. provide a 
nonnegative Nj for all frequencies) only in the low tem- 
perature regime T ^ Tc = \/&P/{7rK), where Rc ^ Rt 
and the iVj-'s of Eq. H54|) are positive quantities for all 
j. On the other hand, in the high temperature regime 
(T > Tc), the solutions Nj provided by Eq. are valid 
only for frequencies ujj ^ tOmax = "H KT{\n£^)/{h y^'), 
where the parameters y^ and ^ are obtained by solving 
(for ^ > 1) the following coupled equations (the first is 
determined by imposing Nj = in Eq. (|54|l . while the 
second is just Eq. (j^ l: 



C^yo - (1 - 7?)^ - = 



j?r(ine) = 



1 — 77 f rjlnS^ 



TT 



Rc 
Rt 



(59) 



with r(a;) defined in Eq. (|49|l . For higher frequencies 
Eq. gives negative values of Nj and we need to choose 
Nj = 0. With these solutions Nj, one can evaluate y^ 
and through Eqs. and ((23) so that from Eq. 
it is possible to obtain 



3 T] ln2 Rt 
R^ 



(60) 



In? 



dx 



9 



- 1 



9 



1 



=2; Tj/y^ 



1 



In the limit T Tc, we find y^ — > ?/ and ^ — > 00, and 
Eq. 1)60(1 reduces to ((57(1 . i.e. the transition between the 
high temperature regime and the low temperature regime 
is continuous. On the other hand, for T = 0, we find 
Rt = 0, so from Eq. ((57(1 we reobtain the lossy channel 
result of Eq. ig^l • The bounds (Ejl and (EOJl are plotted 
in Fig. ^ for different values of T. 



In the low temperature regime we find 

C _ ^ Rt 

^° -[7yi?^ + (l-,)i?|]V2' (55) 
where Rt = J{^^^^ so that the Lagrange multiplier 



VI. DEPHASING CHANNEL 

In this section we will focus on a non-linear noise source 
where the effect of the reservoir depends on the state of 
the message. In particular, we consider the case in which 
the average photon number of the noise source Nj is the 
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Ce 



C 

lower 
bound 



Q 

lower 
bound 



Dj-l 



2(l-'7j)(2Wj+l) 



+ 



1 + 



= e"'^/" f 1 



2(l-,,j)(2Wj+l) 



TABLE III: Lagrange equations for the dephasing channel 
deriving from for the different capacities. 



same as the one (Nj) of the message. The average photon 
number is hence preserved during the transmission. Of 
course, this does not mean that the channel is immune to 
noise: in replacing the lost photons with those from the 
reservoir, some phase correlations are lost. Under these 
conditions, imposing Nj — Nj in @, (|19|l and (|23() . the 
values of the functions ce, k and q become 



ce{N„tj,) = 2[g{N,)-g{{D,-l)/2)] 
HN,,^,) = g{N,)-g{{l-ir,)N,) 
q{N,,Vj) = g{N,) ~ 2g{{D, - l)/2) , 



(61) 
(62) 
(63) 



with Dj = ^1 + 4:N.i{Nj + - tJj), from which the 
Lagrange equations of table IIIII derive to replace those of 
tabled As in the case of the lossy channel the solutions 
of the Lagrange equations depend only on LOj/Q and rj. 
Hence, the same structure as Eqs. apphes, but 

here the functions C, /C and Q are the ones plotted in 
Fig. El 



VII. CONCLUSIONS 

In this paper we extend previous analysis on the ca- 
pacities of broadband bosonic channels with input power 
constraint P, In particular, we analyzed the quan- 
tum capacities in the presence of different noise sources. 
Solutions for the entanglement assisted capacity Ce and 
upper and lower bounds for the classical and quantum 
capacities C and Q were provided. At least in the case 
of unit quantum efficiency (i.e. when the channel is noise 
free), these bounds are tight since they reproduce the 
noiseless capacities Moreover, if the channel noise 

does not have any characteristic frequency (as in the case 
of the loss, white noise and dephasing), the square root 
of the input power dependence (that was known for the 
noiseless case) is reobtained. Even though all the re- 
sults in the paper were obtained by considering a uni- 
form quantum efficiency for all the channel modes, the 




FIG. 6: Plot of the functions C{r]), IC{rj) and Q(r;) for the 
dephasing channel. The hatched region and the cross-hatched 
region are defined as in Fig. Q and contain the classical and 
quantum capacities respectively. The alternative lower bound 
for Q (i.e. C{ri) — 1) is given by the dotted line. As for the 
other channels, it can be shown that Q is null for rj < 1/2 
(dashed line). 



procedure can be extended also to non uniform configu- 
rations. It is also possible to include frequency degerate 
situations, e.g. where one uses polarization degrees of 
freedom to encode information. It is still to be deter- 
mined whether non-linearities in the channel dynamics 
(where some known interaction couples different modes) 
can be used j2| to beat the square root dependence, as in 
the case of the qubit channel discussed in |22(|. 
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APPENDIX A: SINGLE MODE ENTROPIES 

In this appendix we calculate some relevant entropic 
quantities for the single modes when the input is a Gaus- 
sian state. We follow the derivation of Holevo and Werner 
and, for ease of notation, the mode index j is dropped. 

Quantum mutual information. — The quantum mu- 
tual information I{JV, g) for a single mode can be 
evaluated just considering the correlation matrix a of the 
mode input state g, defined as 



(Aq2) i({A(7,Ap}) 
Li({Ap,Ag}) (Ap2) 



(Al) 



where {■, ■} denotes the anticommutator, Aq = q — (q) 
and Ap = p — (p), with q and p the two orthogonal 
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quadratures q = •\/ft/2 {a^ + a) and p = (a^ — a). 

In it has been shown that, for a given value of the 
matrix a, I{JV, g) achieves its maximum value for the 
Gaussian state 



27r 



dz exp 



Aq 

Ap 



(A2) 



where z is a real bidimensional line vector. According 
to Eq. {T)), to determine the value of I{N',g) we need 
the evaluate the input, output and exchange entropies. 
Following ITsL I23 the input entropy of the Gaussian 
state g can be calculated as 



S{g)=g{^X^- 1/2), 



(A3) 



where the function g is defined in (|12() and X± are the 
eigenvalues of a/h. In the same way, we can evaluate 
also the final entropy S{Af[g]). In fact, the state J^[g] 
(evolved by the map J\f defined in (^) is again Gaussian 
and has correlation matrix 



a 



7] a + {1 — r/) B 



(A4) 



where B is the correlation matrix of the jth noise mode 
introduced in (Q. This means that 



S{JV[g])^gi^fy^- 1/2), 



(A5) 



where A'_j_ are the eigenvalues of a' /h. The calculation of 
the entropy of exchange requires to specify a purification 
$g of g: a good choice is the two-mode Gaussian state 



dz / dz exp 



i{Aq, Ap) ■ 



-z(Ag,Ap)-F-(z,z).Af .(2,2)^/2 



(A6) 



where the q and p are quadratures operators acting on an 
ancillary mode and the 4x4 two-mode correlation matrix 
M is 



M 



a (3 
—(3 a 



(A7) 



with /3 = AA/-(A-ia)2 - |-, A being the 2 x 2 matrix 



1 
-1 



The map N ®'t evolves $g into a Gaussian 



state of the same form of (IA6I) with correlation matrix 




-X 



FIG. 7: Plot of the quantum mutual information I[N,q) of 
the single mode as a function of A+ — A- . The different curves 
correspond to different values of the eigenvalues sum A+ + A_ 
increasing from bottom to top. It is evident that the maxima 
are always achieved for A+ = A-. The parameters for these 
plots are iV = .1; r; = .8 . 



where Ai, • • • , A4 are eigenvalues of the matrix A^^M' /K, 

^ A 



A12 being the 4x4 matrix 







In order to evaluate the expressions for the entropies, it 
is convenient to introduce the following real parametriza- 
tion: 



c 



c 

noe" 



(AlO) 



where r is the squeezing parameter. These parameters 
are related through the average number of photons iV by 
the inequalities \/c^ + 1 ^ "0 = (2Af -t- 1 — m)/coshr 
(with m = {q/ft)'^ + {p/H}^): the first relation derives 
from the Heisenberg uncertainty relation, while the sec- 
ond from the energy constraint. The eigenvalues of a and 
a' are respectively 



M' 



a' ^13 



(A8) 



According to 0, IT^ l2^ the entropy of exchange can be 
calculated as 



5((AA® !)[$,]) = i^.g(|A,|- 1/2), 



(A9) 



A± = - [no cosh(r) ± (no sinh(r))2 + c 
A± = ryA± + (l-r/)?i(iV + l/2) , 



■21 , (All) 
(A12) 



while the four eigenvalues of the matrix A-^2^i' /h are 
Ai,2,3,4 = ±[(io ± yi7TT2)/8]i/2, where 
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Lo = 
Li = 



(1 + T]^) - 4(1 - r]fN^ - 4(1 - r7)[l - 7? + i]{X+ + A_)]iV - 2r;(l - 77)(A+ + A_) - 4(1 - TjfX+X^{A13) 
-8(1 - 7/)(l + 2iV)[2(l - r;)(l + 2iV)A+A_ + i]{X+ + A_)] - 4r]\ (A14) 



From these equations and from the definition Q it is 
evident that the entropies of Eqs. (|A3p depend on the 
parameters uq, r and c only through the eigenvalues A± 
of a. Since these quantities are related with the average 
number of photons N by 



A+ + A_ = 2iV + 1 



(A15) 



one can show that the maximum of /(A/", g) for fixed N 
and m is obtained for X+ = A_ (see Fig. |7||. Accord- 
ing to Eq. (|A11|I . this is equivalent to requiring r = 
(no squeezing) and c = (maximally mixed states). Im- 
posing A+ — A- and maximizing with respect to m in 
the above relations, one easily finds that the entropies 
become 



S{Af[g]) 
5((AA(^ !)[$,]) 



9iN) 

giN') 

'D + N- 

9 



N' - 1 



(A16) 
(A17) 

(A18) 



+9 



D-N + N' -I 



where N' and D are defined in (|10|1 and (|11|) . From these 
relations it is immediate to show that the maximum of 
the quantum mutual information for a given value of the 
average photon number N is given by Eq. (j^. 

Coherent information. — Replacing (|A16l) - (|A18p in 
Eq. H21|l allow us to calculate the value of the coherent 
information J{Af, g) for Gaussian states of Eq. ljA2p with 
no squeezing and m = 0. This gives the function q of 
Eq. 123 • 

Holevo quantity.- To calculate the Holevo quantity 
XaA for the code introduced in Sec. ^ we can use the 
above results. In fact both the global state g and its 
components a(^) are Gaussian states of the form l|A2|l . 
with correlation matrices h{N + 1/2)1 and h 1/2 respec- 
tively. Hence it is immediate to calculate the entropies 
that allow to obtain the value of XM reported in Eq. . 
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